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The core of the vortex in the Neel order parameter for an easy-plane antiferromagnet on honey- 
comb lattice is demonstrated to bind two zero-energy states. Remarkably, a single electron occu- 
pying this mid-gap band has its spin fragmented between the two sublattices: Whereas it yields a 
vanishing total magnetization it shows a finite Neel order, orthogonal to the one of the assumed 
background. The requisite easy-plane anisotropy may be introduced by a magnetic field parallel to 
the graphene's layer, for example. The results are relevant for spin-1/2 fermions on graphene's or 
optical honeycomb lattice, in the strongly interacting regime. 



The fact that the low-energy quasiparticles in graphene 
behave as massless Dirac fermions has placed this new in- 
carnation of carbon-based structures at the center stage 
of condensed-matter physics Although in its natu- 
ral state graphene appears to be a gapless semimetal, it 
is tempting to also consider it in its possible strongly 
correlated phases, where the quasiparticles would be- 
come gapped and the ground state insulating. Such a 
semimetal-insulator transition in graphene would in fact 
represent a condensed-matter analog of the Higgs mech- 
anism, with the Higgs field as a quasiparticle composite 
0, 0. Although the ultimate nature of the insulating 
phase that would result from the strong Coulomb interac- 
tion may depend on the details at the scale of the lattice 
constant, the simplified Hubbard model on a honeycomb 
lattice is most likely to have the standard Neel-ordered 
ground state at a sufficiently large U/t. This is certainly 
what is expected at U/t — oo, and is also supported by 
numerical 0] and analytical [3] studies, already for U /t 
above some critical value. 

As long as the quantum phase transition between the 
semimetal and the antifcrromagnctic insulator is not 
strongly discontinuous and not too deep in the insulating 
phase, it is still sensible to consider fermions as itinerant. 
Weak antiferromagnetic ordering is then felt by the for- 
merly gapless quasiparticles as an opening of a "relativis- 
tic" mass gap, with the "mass" being proportional to the 
Neel order parameter. The unusual universality class of 
this transition is completely determined by the presence 
of gapless Dirac fermions on the semimetallic side [3[ . It 
seems natural to ask if there could be other interesting 
consequences of this fact in the insulating phase itself. In 
particular, in light of the known properties of the spec- 
trum of Dirac's Hamiltonian in topologically non-trivial 
backgrounds [f| , [f| one wonders if topological defects in 
the Neel order here may localize sub-gap states at zero 
energy. Such states have been known to display a num- 
ber of unusual features, such as charge fractionalization 
[f| and non-abelian statistics [3] and are being dis- 
cussed intensively in the context of topological quantum 
computing 

I show here that this is indeed the case at least when 
there is an easy-plane anisotropy for the Neel order, when 
the core of the vortex in the order parameter turns out to 



contain a pair of such mutually orthogonal zero-energy 
states. The electric charge of the vortex in this case is 
integer; however, the core states have the positional and 
spin degrees of freedom entangled in an interesting way. 
At half filling, this two-state mid-gap band is half full. 
A single electron occupying the zero-energy band then 
has the direction of its spin alternating between the two 
triangular sublattices of the honeycomb lattice. Its spin 
after averaging over the unit cell therefore vanishes. Its 
Neel order parameter, however, is finite, and points in 
the direction orthogonal to the easy-plane. An electron 
occupying the zero-energy state would thus represent a 
local single-particle antiferromagnet. Placing the system 
in a weak magnetic field parallel to the graphene layer 
is further demonstrated to introduce such an easy-plane 
anisotropy for the Neel ordering. This may help create 
the experimental conditions for an observation of this 
effect in graphene's or optical honeycomb lattice. 

Let us begin by defining the standard Hubbard model 
as H t + Hjj , with 

H t = -t ut(A)v a (A + bi)+H.c, (1) 

H u = Uj2n 1 (X)n l (X). (2) 
x 

The sites A denote one triangular sublattice of the hexag- 
onal lattice, generated by linear combinations of the ba- 
sis vectors a\ — (v3, — l)(a/2), at = (0, a). The second 
sublattice is then at B = A + b, with the vector b be- 
ing either b x = (1/V3, l)(o/2), b 2 = (1/V3, -l)(o/2), or 
63 = (-0/V0, 0). a is the lattice spacing, and U > 0. 

The spectrum of H t for each projection of spin be- 
comes linear in the vicinity of the two non-equivalent 
Fermi points, which we choose to be at ±K, with 
K = (1, l/V3)(27r/av3) [9]. In vicinity of these points 
the linearized Hamiltonian assumes a "relativistically in- 
variant" form Ho = 7o7i<9i with 70 = 1% (8> 03, 71 = 
03 <S> 02, 72 = I2 ® o-\, where {l2,&} is the standard 
Pauli basis. The global "chiral" SU(2) symmetry for 
each spin projection is generated by {735,73,75}, where 
73 = 01 ® 02, 75 = 02 <8 02, and 735 = 17375 = 0-3 ® I 2 
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[lOj. In the continuum limit, the four-component wave- 
function for each spin projection is given by ^>\(x) — 
{u{ a {x),v\ a {x),ul a {^),vl a (x)), where 

Uia{x) = J J^ e ~ lir ' S ^((-y +1 K + 9). ( 3 ) 

and analogously for Vi a (x). The "valley" index i = 1,2 
labels the two Fermi points. The reference frame has 
been rotated so that q x = q ■ K/K and q y = (K x q) x 
K/K 2 : and h = ftB = vp = I, where = is the 

Fermi velocity. 

The interaction term ifj/ can also be written as 

^ = IgElW 1 ) + «(^+ ? )) 2 + «^) - ^(^*+ &)) 2 (4) 

-(/(A) + f(A + b)f (f(A) f(A + b))% 

where n(A) = u\(A)u a (A) and f(A) = 
u\(A)<j at7 ru rJ r{A)i are the particle number and the 
magnetization operators at the site A. Variables at the 
second sublattice are analogously defined in terms of 
v a (B). This is the rotationally invariant version of the 
decomposition previously introduced in Q. 

Assume next that the ground state has a uniform 
average density, zero average magnetization, and finite 
Neel order: (n(X)) = n, (f(A) + f(A + b)) = 0, 
N = (f(A) - f(A + b)) ^ 0. For the Hubbard model 
we believe this to be the leading instability of the semi- 
metallic ground state above critical U/t Q ■ By the usual 
Curie- Weiss decoupling the last term in Eq. (4) yields 
then the mean-field Hamiltonian for the Dirac fermions 
in a fixed Neel background N(x): 

H N =I 2 <E)H - (N(x) ■ a) <g> 7o , (5) 

with the first operator in the direct product acting on the 
spin, and the second on the Dirac indices. For example, 
if N(x) = (0,0,7V) and constant, in the ground state of 
Hjy the average number of spin-up (spin-down) electrons 
is larger on the sublattice B (A). The ground state of 
-ffjv exhibits therefore a uniform Neel ordering along the 
third spin axis in this case. Since all three 7-matrices 
appearing in Hq are diagonal in valley indices, one can 
rewrite Hn = Hi © H 2 , with 

H 1{2) = ±l2®ai{-idi)+h®a2{id 2 )-{N{x)-a)®az (6) 

as the intra-valley Hamiltonians near the Fermi points at 
K (Hi), and -K (H 2 ). Both Hi and H 2 on the other 
hand may be recognized as unitary transformations of 
the generic Hamiltonian H: H = Hq + rh(x) ■ M, where 
M = (i7o73,«7o75,7o)- Specifically, H x = u\HUi, with 

Ux=h® i(?2 (7) 

and m{x) = (Ni(x),N 2 (x),-N 3 (x)), and H 2 = U\HU 2i 
with 

U 2 = ia 2 © I 2 (8) 



and rh(x) = (Ni(x), N 2 (x), ^3(2;)). H represents the 
most general single-particle Hamiltonian in two dimen- 
sions with the relativistic spectrum E 2 = k 2 + m 2 and 
with the chiral symmetry of Ho broken [Til ] . It may also 
be understood as the mean-field Hamiltonian for spin- 
less particles on a honeycomb lattice in the background 
of the order parameter m = (^A/'I'). For instance, 
in = (0, 0, m) is the familiar state with a density im- 
balance between the A and B sublattices [§], 0, [|[, 
whereas rn = (mi,m 2 ,0) represents a state with bro- 
ken translational symmetry, with the "Kekule" pattern 
of hopping integrals of different magnitudes 0] . 

Finding the spectrum of Hm is therefore a special case 
of the general problem of diagonalizing H for a given con- 
figuration of the mass matrix fh(x). Here I focus on the 
issue of zero-energy states for the vortex configuration. 
Assume that one of the components of m(x) vanishes 
everywhere, say m^(x) = 0. As a result, {70, H} = 0. 
Since Jq = 1, 70 = P+ — P- where P± are the projec- 
tors onto two orthogonal eigenspaces corresponding to ±1 
eigenvalues. This, on the other hand, means that when 
m 3 (f) = 0, H = P+HP- + P-HP+, so that H is block- 
off-diagonal in the eigenbasis of 70 in which P + = I 2 © 
[i^ . For 70 as defined right below Eq. (1) this change 
of basis amounts to a simple exchange of the second 
and the third components of ^>(x). Next, assume that 
A = nii(x) + im 2 (x) = |A(r)|e lp *, p = ±1, is the vor- 
tex configuration, with (r, <p) as polar coordinates and, 
and |A(r — > 00) | — > m. This problem was considered 
by Jackiw and Rossi [l4| and recently by Hou et al 0]. 
Defining d z = d x + id y , the zero-energy state ^o{x) sat- 
isfies either 

id z v(x) + A*v*(x) =0, (9) 

with vi = v 2 — v(x) and ui = u 2 = 0, or 

id g u(x) + A*u*(x) = (10) 

with ui = u 2 = u(x), and vi — v 2 — 0. For p = 1 one 
finds: 

v(x) = -e-^fir), (11) 
r 

with u = and - ln/(r) = ^ \ A(t)\dt + i(n/4) , and 

u(x) = Cf(r) (12) 

with v = 0. C is the normalization factor. Only the 
second solution is normalizable, however, and thus rep- 
resents the unique zero-energy eigenstate of H . The sit- 
uation is reversed for an antivortex, when p = — 1. For 
general p one can show that there are \p\ linearly inde- 
pendent zero-energy bound states. 

An important general property of the zero-energy 
states should be noted. Consider the following sum: 

q(x) = ^ e (S)Q^e(x), (13) 
EeR 
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where Q is a traceless Hermitean matrix, and {^e(x)} 
the eigenstates of the generic Hamiltonian H . If the sum- 
mation is performed over the whole spectrum, q(x) = 0. 
If R includes only the occupied states, on the other hand, 
the sum represents the ground-state average of a physical 
observable, (q(x)). By subtracting a half of the vanishing 
sum over the whole spectrum one can rewrite this same 
average as [lj| 



<«(*)> = 5(£- E 

occup unoccup 



)*yf)Q* B (f). (14) 



The last form makes it evident that if there exists a uni- 
tary matrix T that anticommutes with H, for any Q that 
commutes with T the only states that can contribute to 
(q{x)} are the states with zero energy. For the vortex 
Hamiltonian considered above, T = 70. Choosing Q = 70 
yields then (q(x)) = (m 3 (x)), whereas for Q = 77073, 
(q(x)) = (mi(x)} for example. The important observa- 
tion is that the contribution to the hard-axis component 
of the order parameter comes exclusively from the zero- 
energy state. 

We now translate these results to the vortex in the 
Neel order. Take the Neel vector to be in the plane 1-2: 
N = (jyi(f),JVa(5),0), with A(x) = N^x) + iN 2 (x) = 
|A(r)|e^. The zero-energy eigenvalue of H\ is then 



#1.0 = 




= f(r) 



(15) 



and of H2, 



*2,0 = 




= U^ 



f(r) 



/o 

i 

Vo 



(16) 



with f(r) assumed to be normalized. The two zero- 
energy states are degenerate, orthogonal, and exponen- 
tially localized at the center of the vortex. Consider 
a linear combination \jd = a*(*| ,0) + &*(0, ^ ), 



\b\ 2 



1. In this state, 



u T (x) 



ae 



iK-x 



and 



(17) 
(18) 

(19) 
(20) 

v. Introduc- 
ing tv — (ut,u*A and — (vtjVV), the average particle 



H (x)=ibe- iK - s! f(r), 

vi(x) = be- iR - £ f{r), 
v l (x)=iae i *-*f(r). 



where w a — wi a e tK ' x + W2a& x , w = u. 



densities on the sublattices A and B in the state \P are 
equal to 



u)u = v'v = |/(t") I"" 



(21) 



independently of the state \F The average of the third 
component of the spin, on the other hand, is 



u=-v^ 3 v=(\a\ 2 -\b\ 2 )\f(r)\ 2 , (22) 



u'a 3 



and thus alternating between the two sublattices. For 
the easy-plane components of the average spin we can 
compactly write, 



iu'a+u = iv^<T-v 



2a*be 



-2iKx 



(23) 



which is evidently zero for both ^ = 4"i and \I/ = \E f 2 

Exactly at half-filling the two-state band at zero- 
energy is occupied by a single electron, which is in some 
arbitrary state Whereas in any such state the elec- 
tron has an equal probability to be found on either sub- 
lattice, the projection of its spin along any axis when 
integrated over the whole system is zero [131 ■ However, 
unless \a\ — \b\, the electron spin will manifest itself in 
the staggered magnetization in the direction orthogonal 
to the plane of the Neel background. Taking into account 
the zero-states' local contribution to the order parameter 
the vortex would thus be turned into a half-skyrmion. It 
is a remarkable property of the zero-energy states that 
the spin of a single electron is distributed in such an al- 
ternating pattern. It agrees, however, with the observa- 
tion on the general nature of the zero-energy states made 
earlier. 

As shown below, one way of introducing such an easy 
plane for the Neel order parameter is to turn on a uni- 
form magnetic field and couple it to particle's spin. It is 
of interest therefore to understand how that would per- 
turb the zero-energy states. In the zero-energy subspace 
spanned by (\1/} ,0) and (0,^20) the Zeeman term is 
represented by the perturbation 



Hz = A((T 3 ® I 2 ) 8 (er 3 ® h), 



(24) 



where A = —(ibB. It is obvious that ("J^gl-Hal^o) = 0, 
and there is no shift in energies to the first order in B. 
Physical reason is precisely the states' vanishing magne- 
tization, which decouples them from the magnetic field 
to the first order. 

The energies do shift, however, if the Neel vector is 
tilted out of the easy plane. Assume a configuration N = 
(JVi(f), N 3 (a!),Ns), with N 1 (x)+iN 2 (x) = |A(r)|e* as 
before, but with a uniform N 3 -C |A(oo)|. Within the 
two-dimensional zero-energy subspace this perturbation 
is represented by 



-N 3 (a 3 (g>a 3 )®{a 3 (g>a 3 ). 



(25) 



This gives (vPj _Q\H 3 \^j ) — ±A 3 (5.y. A finite component 
A3 at half filling would therefore force one of the states 
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in Eqs. (15)-(16) to be occupied and the other one to be 
empty. 

Let us show how the Zeeman coupling to the magnetic 
field introduces easy-plane anisotropy for the Neel order 
in the present case. To this purpose add the term A(o3 <g> 
I) to Hn in Eq. (5), and assume a uniform N. The third 
spin axis, of course, may be in an arbitrary direction in 
the real space. Decomposing the new Hamiltonian again 
as right above Eq. (6), one now finds Hi and H-> as 
unitary transforms of H + A735 , with the eigenstates [18| : 

E(k) = ±[{^Nl + k^± |A|) 2 + (N 2 + Nl)} 1 / 2 . (26) 

In the ground state at half filling all the negative energy 
states are occupied and the positive energy states empty. 
It is easy to show then that among all possible orien- 
tations of N the lowest energy belongs to the one with 
-/V3 = 0. An obvious way to enforce such a purely Zeeman 
coupling would be to orient the magnetic field parallel to 
the plane of the honeycomb lattice. The Neel order pa- 
rameter would then be confined to the plane orthogonal 
to the honeycomb lattice. 

Hou et al. Q have recently found zero-energy states on 
graphene's honeycomb lattice in the spectrum of spinless 
electrons in the vortex in their two-component "Kekule" 
order parameter, as discussed below Eq. (8). In that case 
there is a single zero-energy state, a remarkable conse- 
quence of which is that the charge of a vortex that binds 
it is fractionalized. In contrast, since we have two zero- 



energy states, the charge of each is simply unity, but it is 
the spin properties of the states that are nontrivial. The 
crucial feature of our example, however, is that the rele- 
vant U(l) symmetry of the order parameter is the exact 
rotational symmetry. This makes one cautiously opti- 
mistic about the survival of the mid-gap states once the 
lattice effects are fully restored. In this context it would 
be particularly interesting to diagonalize the lattice ver- 
sion of Hn numerically. A related recent calculation [l9j 
shows zero-energy states to be surprisingly resilient to 
the effects of the discrete lattice. 

In conclusion, the mean-field Hamiltonian for Dirac 
quasiparticles in the background of a vortex in weak Neel 
order on graphene's honeycomb lattice has two orthogo- 
nal core states at zero energy, with positional and spin 
degrees of freedom maximally entangled. An electron in 
one of these states at half filling is a single-particle an- 
tiferromagnet with a finite, localized contribution to the 
staggered magnetization in the direction orthogonal to 
the easy plane, and zero average magnetization. The re- 
quired easy-plane anisotropy was shown to be introduced 
by placing the graphene layer in a weak parallel magnetic 
field. 
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